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1. Introduction 

Supersymmetry is a beautiful and phenomenologically interesting candidate for phys- 
ics beyond the standard model. Our understanding of quantum field theory has been 
improved by our ability to make exact statements within supersymmetric theories, 
such as non-renormalization theorems and dualities. 

There is a large class of nonperturbative results or conjectured results for super- 
symmetric theories, particularly supersymmetric Yang-Mills theories, both in four 
dimensions []H ^|, |3| and in three dimensions 0, [5|, H |7j. For instance, in 4 dimensional 
N=2 super Yang-Mills (SYM) theory, Seiberg and Witten's work || essentially solves 
the theory. In three dimensions, a number of exquisite results have been obtained 
involving mirror symmetry in N=4 SYM theory |3J] and also in N=2 theories §]. 
In [|7| these symmetries were motivated in the context of string theory. There is also 
a beautiful conjecture by Witten M which makes N=l SYM theory in 3 dimensions 
an excellent testbed for spontaneous (dynamical) supersymmetry breaking. 

It would be very helpful if we were able to test some of the techniques for studying 
supersymmetric theories, for instance by "solving" the theories involved in a non- 
perturbative way. The lattice is the best candidate method for solving such theories. 
Unfortunately, the lattice regulator almost inevitably breaks the supersymmetry. 
The lattice is, after all, a regularization scheme designed to preserve exact gauge 
symmetry at the expense of manifest Poincare invariance. Since supersymmetry is 
a space-time symmetry, i.e. an extension of the Poincare group, it is of no surprise 
that it is broken on the lattice. 
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It is certainly possible to use the lattice to study a theory which possesses a 
symmetry broken by the lattice action. After all, the main application of the lat- 
tice technique is the study of QCD, which has 0(4) Euclidean invariance, broken to 
(hyper)cubic symmetry on the lattice. In this case, the symmetry is recovered auto- 
matically, up to corrections suppressed by the lattice spacing, because the low-energy 
phenomena are described by an effective theory which possesses full 0(4) symmetry 
accidentally. That is, the low dimension operators which respect (hyper)cubic sym- 
metry also respect 0(4) symmetry, so that symmetry is recovered in the infrared. 
Unfortunately, for supersymmetric theories containing scalars, the scalar mass term 
is low dimensional and breaks the supersymmetry. Therefore very fine tuning of the 
lattice action is generally required to recover supersymmetry in the infrared. 

A great deal of work has gone into looking for a way to preserve supersymme- 



try, at least in the small lattice spacing limit |], [II], [12], [II], [L3|]. In 4 dimensions, 
it appears that one should be able to implement N=l SYM theory without matter 
[[14|, [L5|, jlS] , that is, a theory containing only gauge fields and an adjoint representa- 
tion, Weyl fermion (the gluino). This is because this case does not involve a scalar 
field, so supersymmetry is an accidental symmetry once enough chiral symmetry is 
present to prevent a mass for the gluino. This is now possible with recent advances 



in 4-D chiral fermions |17, 18, 19, [20||. For cases involving scalar fields, most at- 



tempts involve building in enough supersymmetry into the lattice construction to 
ensure that scalar mass renormalization does not occur. This is true, for instance, 
of proposed implementations of the Wess-Zumino model in 2 dimensions ||, H)), ^1 



Various 2 and 3 dimensional, highly supersymmetric pure Yang-Mills theories can be 
implemented by a technique developed by Kaplan, Katz, and Unsal []T32[ (see how- 
|23|). A technique for putting N=2 SYM theory on the lattice using "twisted" 



ever 



supersymmetries with Kahler-Dirac fermions has been developed by Catterall p4 



and recently extended to N=A supersymmetry in 4 dimensions |25fl . Also D'Adda 
and collaborators [26 have recently completed the construction of a 2-d lattice which 
retains all the supersymmetries of the target SYM theory exactly on the lattice us- 
ing twisted supersymmetry and "fermionic" links where supercharges live in direct 
analogy to gauge fields living on bosonic links. This would seem to be well motivated 
by analogy to the continuous superspace formulation of supersymmetry in which the 
generators of the supersymmetry are regarded as translation operators along new 
Grassman valued coordinates. 

These techniques often come with certain limitations such as confinement to 
even dimensions, to very highly supersymmetric theories or to specific numbers of 
fermionic fields. They are frequently made difficult by highly non-trivial lattice 
actions that can contain a wide range of complications such as non-commutativity, 
the failure of reflection positivity, unphysical moduli, or sign problems with complex 
fermion determinants. Most of these issues can or have been overcome to varying 
extents, yet they still remain to complicate the treatment. 
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We argue here that, for all the 3-dimensional theories of interest with N=2 (4 real 
supercharges) or more, one can proceed by actually doing the fine tuning on the most 
conventional lattice action available, the Wilson action with Wilson fermions. The 
reason this approach is available is because these theories are super-renormalizable. 
That means that loop corrections involving the ultraviolet converge in powers of the 
lattice spacing, so the tuning only requires a lattice perturbation theory calculation 
to a finite loop order, which turns out to be 2 loops. Further, we explicitly do this 
calculation for the case of SYM theory plus fundamental matter. This holds out the 
possibility of testing some very interesting claims || |J of exact results in N=2 SYM 
theory with matter. Some similarly spirited work in the context of 2-D N=2 pure 
SYM has recently appeared in . 

In Section 0, we will explain the general idea which makes it possible to preserve 
supersymmetry, in the small lattice spacing limit, in 3-dimensional, N=2 supersym- 
metric theories. In Section [I] we present in detail the implementation of the Wess- 
Zumino model, that is, a theory of fermions and scalars. In Section |] we present 
the implementation of gauge theories. Finally, some concluding remarks appear in 
Section |5|. 



2. Super-renormalizability and SUSY breaking 

No one is interested in the results of a lattice calculation per se. After all, a lattice 
"field theory" is actually just a statistical model, not a true field theory. The reason 
that the lattice technique can teach us something about field theory, is that in the 
infrared, the correct effective description of a lattice theory IS clS cL quantum field 
theory. What one must do is to ensure that the infrared behavior of the lattice 
theory coincides with the continuum quantum field theory of interest, or at least 
that it does so in the small lattice spacing limit, so that the behavior of the field 
theory can be probed by making a zero lattice spacing extrapolation. 

It is easy to write down a lattice gauge theory which, at tree level, will look in the 
infrared like the theory of interest 1 . The problem is that, in the UV (at the lattice 
spacing scale), the lattice theory typically does not have the full symmetries of the 
theory we are interested in. Generally it is possible to formulate lattice theories so 
that they have exact gauge and (hyper)cubic symmetries. However, under super- 
symmetry, the variation of a fermionic field can involve the derivative of a bosonic 

1 In four dimensions this statement is true of vectorlike theories, but serious complications arise 
if one wants a chiral theory, that is, a theory with two component spinors in a representation which 
is not real or pseudoreal and which are not balanced by an equal number of spinors of opposite 
handedness in the same representation. To date, it is understood how, in principle, to put certain 
Abelian chiral theories on the lattice; it it not clear, even in principle, how to put more general 
theories on the lattice while still retaining manifest gauge invariance |^8| . Fortunately, for the three 
dimensional theories of interest to us, this is not relevant, because chiral fermions do not exist in 
odd dimensions. 



3 



field; and since derivatives become finite differences on the lattice, supersymmetry 
will generically be badly broken at the lattice spacing scale. Furthermore, even if 
we construct the lattice theory to satisfy supersymmetric relations in the infrared, 
radiative effects involving UV (SUSY breaking) modes will typically communicate 
those effects to the infrared modes of interest. 

The IR effective theory is not the tree level theory. Rather, it is the theory 
one obtains, by writing down the most general continuum quantum field theory 
consistent with the field content and symmetries of the lattice, and performing a 
matching calculation between the lattice theory and that continuum effective theory, 
to determine what the actual parameters of the IR effective theory are. For instance, 
if we made a tree level lattice implementation of the Wess-Zumino model (which 
in 3 dimensions exhibits N=2 supersymmetry, that is, it has 4 real supersymmetry 
generators) , 



with $ a complex scalar and ip a two component spinor, then we would generically 
recover an infrared theory where all terms permissible with this field content were 
present; 



Here and represent the difference in field normalization between the lattice 
and continuum fields; they can be removed by a field rescaling, but we must keep 
them in mind when we compare lattice correlation functions with their continuum 
counterparts. The point is that the IR behavior typically involves radiatively gen- 
erated terms which do not respect the intended supersymmetry. In particular one 
does not expect m 2 ^ = m 2 . 

In 4 dimensions this problem is severe. The SUSY violating, radiatively induced 
terms appear at all orders in perturbation theory, with coefficients, at high order, 
which are only suppressed with respect to the lower order coefficients by powers of 
a dimensionless coupling. Further, additive scalar mass renormalizations are diver- 
gently large at every loop order. That is, in 4 dimensions, the contributions to the 
mass squared parameter are of order 



5m 2 at 1 loop: \ 2 /a 2 ; 2 loops: \ 4 /a 2 ; 3 loops: \ 6 /a 2 ; . . . , (2.3) 



where a is the lattice spacing. Every such coefficient is problematic; a severe nonper- 
turbative tuning is needed to remove them. It is not at all clear how to perform such 
a tuning; generally we can only perform nonperturbative tunings in lattice gauge 
theories if we have one exact conservation law or Ward identity per tuning required. 




(2.1) 



C m = Z^*d^ + Z^frl) + mj$*$ + m^V 

+ (\ y ®il> T ei/> + h.c) + \ 2 s U*<f>) 2 + (High Dim.) 



(2.2) 
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The beauty of 3-D is that the desired theory is generally super-renormalizable. 
Consequently, the UV is very weakly coupled; specifically as the lattice spacing is 
taken to zero, the coupling at the scale of the lattice spacing falls linearly with lattice 
spacing a. This means that, while the SUSY breaking nature of the UV regulator 
radiatively induces SUSY breaking effects in the IR, the matching calculation which 
determines them converges very quickly. At each loop order, we determine the match- 
ing of parameters to one more power of the lattice spacing a. For instance, in the 
above model, if we compute the mass squared for the scalar field, generated by UV 
physics, the contributions at different orders in the loopwise expansion are again of 
order A 2 , A 4 , A 6 , .... But A 2 has mass dimension 1. Since the matching calculation 
involves only UV physics, the only scale which can balance the explicit powers of 
mass is the lattice spacing scale. Therefore, the terms in the loopwise expansion are 
of order 

8m 2 at 1 loop: A 2 /a; 2 loops: A 4 ; 3 loops: aA 6 ; .... (2.4) 

The one and two loop contributions are significant and must be removed by an 
appropriate counterterm. However, three and higher loop effects vanish in the a — > 
limit, and so can be neglected. For the scalar self-coupling A 2 , the one loop correction 
is already 0(aA 4 ), and so a tree level treatment is already sufficient. 

We should stress here that the infrared cancellation which ensures that only 
the scale a can appear to balance the explicit mass dimensions in the coupling, is 
a generic property of matching calculations in effective field theory and contains 
no statement about the IR physics of supersymmetry. Infrared divergences arise 
from large length scale (or low loop momentum) behavior. By construction the two 
theories for which any matching in an effective field theory formulation is being 
performed have the same IR behavior, and so any IR divergence will cancel in the 
difference. The statement is then simply that the IR extension of the lattice theory 
does describe the SUSY theory of interest provided that it contains the same degrees 
of freedom and that the coefficients of the terms in the Lagrangian can be matched 
with those of the theory of interest. In the 4-d case this matching does not seem 
possible within a perturbative framework; however, in a super-renormalizable theory 
it certainly is. 

We see that only a finite loop order is needed before all remaining corrections 
are suppressed by powers of a. It is therefore feasible to perform the matching 
calculation to the requisite order analytically, and to tune the lattice theory based 
on the purely analytic result of this perturbative matching calculation, to ensure that 
the IR effective theory satisfies all relations implied by SUSY up to a suppressed 
corrections. 

How hard is this tuning? In the present paper we will be satisfied with removing 
SUSY violating effects which do not vanish with the lattice spacing; that is, we 
will leave 0(a) SUSY violating effects, but prevent 0(a°) SUSY violation. In this 
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case it is straightforward to show that the scalar masses must be determined at 
two loops, and the fermionic masses must be determined at one loop. All other 
Lagrangian parameters may be taken at tree level. This same technique has been 
applied to purely bosonic 3-dimensional Yang-Mills Higgs theory, in the context of 
understanding the electroweak phase transition, in |29|, |3(| |31], |32], Note that 
we will not attempt to provide renormalized operators, and in particular we will not 
compute the SUSY violating renormalization of the vacuum energy, which would 
require a 4 loop effort. 

To conduct the matching calculation, we need to compute some set of correlation 
functions in both the lattice theory and the continuum, supersymmetric theory, and 
equate the answers. At the level of interest, we need only do so for two correlation 
functions, which must be sensitive to the scalar and fermionic masses. The obvious 
candidates are the respective two-point functions (self-energies). Further, there is no 
need to compute the supersymmetric values, since we already know that they vanish 
identically. Therefore, all that is required is to compute the fermionic self-energy 
at zero momentum, at one loop, and the scalar self-energy at zero momentum, at 
two loops, and to assign counterterms to cancel these contributions. Provided that 
we take suitable combinations of loop contributions (fermionic and bosonic), we are 
guaranteed that all loop integrations will be IR finite-precisely because supersym- 
metry is satisfied in the infrared. Therefore, the matching calculation will consist of 
computing a handful of IR finite linear combinations of Feynman graphs involving 
scalars, fermions, and gauge fields. 

In our calculations here, we will eliminate all corrections which are unsuppressed 
by powers of the lattice spacing a; but corrections suppressed by a single power of a 
will still exist. If we were very strong, it would be possible to eliminate 0(a) errors 
as well, by conducting a 1 loop determination of couplings, a 2 loop determination 
of fermionic masses, and a 3 loop determination of scalar masses. (Tree level im- 
provement of the fermionic action would also be required.) Such a calculation has 



been carried out for three dimensional, O(N) symmetric scalar field theories [[34 



indeed, in this case many of the 0(a 2 ) corrections have also been eliminated ||35|| . In 
principle it is possible to eliminate errors to any fixed order in a, by using a suitably 
tree-level improved action and evaluating a number of low dimension operators to 
suitable loop levels; as the number of powers of a accuracy desired increases, more 
and more operators must be improved. We will not attempt such a program here. 



3. Detailed treatment of the Wess-Zumino model 

Consider a theory of scalars and fermions with N=2 supersymmetry in 3 dimensions. 
The matter content is the same as for a 4 dimensional N=l supersymmetric theory. 
To make the extension to gauge theories simpler, we consider a theory with a global 
SU(iV c ) symmetry, and three sets of fields; those in the fundamental representation, 
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$/, ipf] those in the antifundamental representation, $ a , ip a ; and those in the singlet 
representation, $ s , ip s . 

The most general superpotential, excluding mass terms, is 



W — A u<& f $ $ !. + — $ $ $ i, 



(3.1) 



with C,ijk totally symmetric in its indices. We have left out super symmetric mass 
terms simply because including them will not lead to any change in the mass coun- 
terterms we will need; any loopwise correction involving a mass term will be 0(m 2 ); 
at worst this leads to 0(m 2 \ 2 a) effects, which are beyond the order under consider- 
ation. 

The Lagrangian which follows from this superpotential is 



(d^f^id^fj + 4, ftfo + (/ - a, S)] + \* jk \ilmK,f ®a,l$:, k $s, m 
+\* ik \ lirn $l j - <S> f;l <S>* >k $ s , m + \) ki \ lmi $>l r K,k®f,r ®a,m 

<I> /./ <1> "./ <1, :./ <1> ; 



+ 



+ A 



Hjk 



+K 3k \*wh e + *MU e raj + KA, k e r u 



(3.2) 



where the dot product • refers to SU(iV c ) indices and e is the antisymmetric 2x2 
matrix, e = ia 2 . Note that the 3-dimensional gamma matrices are 2x2 matrices, 
and are in fact nothing but the Pauli matrices, a^. 

The lattice implementation of this theory is obtained by discretizing space onto 
a cubic lattice of spacing a, replacing f d 3 xC with a 3 J2 X C, and writing the gradient 
terms as follows: 



/ 



d 3 x (d fl $y(d»$) 



d 3 x ^ 



x,n 

a 3 ^ 



§*(x+a£i) - <&*(z) $(x+ap,) - $(x) 



(3.3) 



X,fl 



^(x+ajj) -j)(x-aji) 
2a 

ar —ip^x+afi) + 2ip(x) — ip(x—afi) 



• (3-4) 



The "extra" term in the fermion gradient term, which for slowly varying fields looks 
like — (ra/2)^d 2 ^, is called the Wilson term, and is required to remove fermion 
doublers. It vanishes in the continuum a — > limit, but its presence complicates 
the treatment, particularly because the value of the Wilson coefficient r is to be 
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chosen by the practitioner (though r > is required and r < 1 is desirable to ensure 
reflection positivity). We will present results for the two values of r = 1 and r = |. 

In addition to these, it is necessary to add the following mass (counter) terms to 
the lattice Lagrangian: 

+6M Uj ^ fd + SM^lrtoj + SM.^lrf.j . (3.5) 

In principle we should also allow for multiplicative renormalization of the couplings 
and wave functions by including Z factors for each, but this is not needed at the 
order of interest. 

To write the propagators of the lattice fields it is convenient to define the quan- 
tities 



Z^~ sm 9 ' P l ^ = -smp fM a, (3.6) 



p 2 

■ a 2 "~ 2 ' ' " a 



in terms of which the propagators of the scalar and fermionic fields are 

A p = ^> S p = 2 ■ ( 3 - 7 ) 

p l p l + M 2 

Here M p = ^-p 2 is the momentum dependent effective lattice mass induced by the 
Wilson term. 

We need to do three integrals using these lattice propagators to complete the 
matching calculation; define 



C ys _ f d 3 p 



47ra ~ J BZ (2tt) 3 

Cyf __ f (1 3 P 

4vr " J BZ (2tt)3 

167T 2 



{2A p +TiS p S p y (3.8) 
A P (S P + S^ P ), (3.9) 



J {^k\\+k - A fc Tr (Sp(S k +p - S p ) + Sl +p (S p - S k+P )^ 



(3.10) 

Here BZ means that the integral is to be taken over a cubic Brillouin zone of extent 
(— 7r/a, 7r/a] 3 , that is, each momentum component runs from —ir/a to n/a. We have 
included the factors of Air to imitate the expected behavior that each loop order 
involves y 2 /47r, and to match notation with previous literature [29|. The diagrams 
responsible for the contributions are presented in Figure [|; the integrations for C yy , 
the two loop scalar mass correction, have been written in the way presented so that 
the integrand is absolutely integrable and carries all of its singularities on the three 
surfaces k = 0, p = 0, and p+k = 0, which is convenient for numerical evaluation. 



S 




Figure 1: All diagrams needed in scalar and fermionic mass corrections; A are for C ys , B 
is for Cyf, and C are for C yy . Thin lines are scalar s, heavy lines are fermions. All of the two 
loop diagrams have the same combinatorial factor, so the integrands can be added directly 
to produce an IR finite overall integrand. The heavy dots are one loop mass counterterms, 
given by the one loop diagrams evaluated at vanishing external momentum. 



We can readily verify that, at low momenta, each integrand is well behaved. At 
leading order for small p, the integrand for C ys behaves as 

* Tr (ip 1 ) 2 



2A(p)+ Tr S 2 (p) ~ — + 



p" 



2 2 

= — - — = (each term is +0(p°)) . (3.11) 

p2 p2 



Similarly, the integrand of C y f goes as ra/p 2 at small p. The integrand for C yy 
requires more patience to analyze, but it also proves to be well behaved for p — > 0, 
for k — > 0, and for p, k simultaneously taken to zero. 

None of the integrals can be done in closed form (to our knowledge), but all are 
relatively tractable by quadratures; we find numerically, 

C ys (r = 1) = 6.4706034146527591308 C ys (r = 0.5) = 5.057247581039541 
C yf (r = 1) = 2.29977456857632 C yf (r = 0.5) = 2.22804716126902 (3.12) 

C yy (r = 1) = 5.425954134(5) C yy (r = 0.5) = 6.8513618(8) . 

In terms of these coefficients, the required renormalizations of the masses are 



<WU(! loop) 



a 



ys 



Area ' 



Sm 2 Uj (l loop) = -\i lm \jlr 



a 



ys 



Attq ' 
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lo °P) = -Km^ljmj^- , (3-13) 



for the scalar masses at one loop, 



5M 8)i j — yd { X* lmi Xi mj + ~£*im€jlmj > 
Q 

JM^ = A* m A ym ^ , (3.14) 



for the fermions, and 



^ m s,y(2 loop) — | + diKimi^nqjX*i qk \l m k + dfKnni^nj^q'lk^rnlk 

df£, qn i£,qmj ^klm^kln + ^kmSlrnj^knq^lnq | 
5m /,ij( 2 l0 °P) = { A *mn\gnAfc^A fcm/ 



16tt 2 



1 l c. 

dfXi nrn Xjnq^klq^klm + Ti^ikl^jkra^nqm^nql 



J yy 



^klq^klm T ^ 'Hkl'^jkm^nqm^nql J ^g_2 ' 
£wii( 2 loop) = {A^A^nA^Amfcj 

1 ~\ c 

+ di^nim^njq^lkq^lkm + g A&« ^kjm£ nqm £nql J ' (3-15) 



for the scalar masses at two loops. Here df is the dimension of the fundamental 
representation, d{ = N c in SU(iV c ) gauge theory. The full scalar mass counterterm is 
the sum of the 1 and 2 loop contributions. 

This completes the renormalization of the theory at a level which will leave only 
0(a) supersymmetry breaking effects. 



4. N=2 SU(iV c ) gauge theories with fundamental matter 

Now we extend these results to the case where gauge interactions are also present. We 
will not review here, how gauge fields are put on the lattice; we refer the interested 
reader to H. Rothe's book f3~6fl , which also presents the Feynman rules associated 
with the pure gauge and gauge/fermion sectors with 1/4 — ► C^/12 in Eq. (15.39) 
and 2/3(5 A b5cd + ...)-► J2^(T A T B T C T D ) in Eq. (15.53b) for the generalization 
to SU(iV c ) (the sum is over all permutations of ABCD). 

The added fields are a gauge field A^, an adjoint fermionic gaugino %■> an d an 
adjoint scalar field we will write (p. In the approach where we obtain a 3-D N=2 
super symmetric theory by dimensional reduction of a 4-D, N—l SUSY theory, the 
scalar is the gauge field component in the direction which was compactified. Besides 



10 



their gauge interactions, which follow from each species' gauge representation, these 
fields also introduce new Yukawa and scalar couplings, 

£(new) = V2g^\.e(x A TT A ^ u + ^ u T A ^J.ex A + (f^a)) 

+ 9 j($%T A <S> u + SVT^) ($*fjT A $fj + KjT^a,) 
+g 2 <P A <P B ($* Ll T A T B <& } , + $* ii T A T B $„ |i j 

+g<f{f Jt F A il> J ,i + <,T A ^) + #V^X, (4.1) 

where = —iJabc is the generator of SU(iV c ) in the adjoint representation and T A 
must be treated in context as the group generator in the fundamental representation 
when between / species, or the antifundamental representation (— T*) when between 
a species. These interactions, which are local, can be implemented on the lattice 
in the obvious way. The gauge interaction implementation is not so simple but has 
been well treated in the literature. 

Though the matching calculation as a whole is perfectly IR safe and gauge invari- 
ant, individual diagrams are, in general, not. The calculation reduces to determining 
IR safe combinations of diagrams that are gauge invariant. To this end we define the 
invariants of the group representation as 

Tr (T A T B ) = T F 5 AB 
(T A T A ) ab ee C F 5 ab 
Tr (F A F B ) = {F C F C ) AB ee C a 5 ab . (4.2) 

Here Tp=C (N c ) is the first invariant of the fundamental (defining) representation of 
SU(iV c ) also called the trace normalization. It is usually chosen to be \. With that 
choice the second invariant is C F =C 2 (N c ) = (N 2 -l) /2N C . C A =C(G)=C 2 (G)=N C are 
the invariants of the adjoint representation of the group, rif and n a are the number of 
fundamental and anti-fundamental multiplets respectively. Each diagram will include 
a prefactor involving some combination of these invariants, the Yukawa couplings, 
the dimension of the fundamental representation and the sum rif + n a . Since the 
matching calculation is valid for any Lie group, we can factorize the diagrams based 
on this prefactor. 

Two other group theoretic relations are needed to complete the calculation, which 
we provide here since they appear in the results almost unchanged: 

rpBrpArpB = ^ _ JlA ^ F^Brj^C = ^C A T A . (4.3) 

New counterterms are needed, which we name \5m 2 ^ and SM X ; they have the 
obvious definitions. The new contributions to the one loop mass counterterms turn 
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out to require only one new lattice integral, 



^W_^ + I^!/^,J (4 . 4) 



4tt J bz (2tt) 3 V 2 p 4 \p 2 + M 2 
which is known analytically for r = 1 and easily determined for r = ~: 

C ff/ (r = 1) = -S/2 C 9/ (r = 0.5) = .097938749331668 (4.5) 
where £ = 3.17591153562522 [||]. They are, 



Ana 



( 5m 5).4i? = -9 2 {i n S + n a )T F + C A y 
(5M aJ ) ab = g 2 C F 8 ab ^-, 



Cys 

Ana 



(5M X ) AB = g 2 C A ^- + g 2 ((n, + n a )T F + C A )«5^^ . (4.6) 

The new two loop counterterms come in many combinations factorized according 
to the prefactors determined by the group theory (i.e. products of the group invari- 
ants). In terms of the nine new integrals defined in appendix [A] the new counterterms 
are (obviously X* lm Xji m — > X* im Xij m for 8m 2 . in the second expression) 

Smlij = g 2 Kmi^ m jC F dF^r^ , (4.7) 

Cfund ( y^ifund 

„1 a - - I _ - - . < ~~'g2 



i Sm lf ;ij) ab = 9 2 8abKlm X jlmC F -^r^- + g%jSab\ T F C F (n f + U a 



167T 

' C fund E2 \ / C /W 2 



2 



16?r 2 3 16vr 2 j I 16vr 2 18 16tt 2 

1 »' I - - • ( - - U g3 

167T 



2 



+ (C A ) 2 I - — ^ 1 - -T F C A ( C F - -C A ) ^ • (4.9) 

1 16tt 2 1816vr 2 J 3 V 6 / 16tt 2 J v ; 

The last term in each scalar mass correction is from the piece of the 4-point gluon 
vertex with a group structure that is fully symmetric between the four lines. It is 
separately gauge invariant and IR finite, and (as shown) can be found in closed form 
in terms of the constant E, defined in Eq. |A.6j . The group theoretic factor preceding 
it arises as 8 AB 8 CD E T A T B T C T D . 



C fi 5 S 2 \ 4^ / 1 W47TE 
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The constants appearing in these expressions are 



C } g T\r = 1) = -2.647013(1) 
C g f j (r = 1) = 7.75588650(2) 
C g f(r = 1) = 17.536258926(7) 
C a J j (r = 1) = -.3347923(2) 
C°*'( r = 1) = 13.0938429(1) 



C s g in9 (r = 1) = 3.588328893(6) 
C^ ,nd (r = 1) = -4.89236097(1) 
C{l nd {r = 1) = 10.2296763(2) 



C'5 ,nd (r = 1) = 22.712647140(8) 





(4.10) 



C"f (r = 0.5) = 20.658655(8) . 



The number in parentheses is the error in the last digit. It was determined conser- 
vatively since the accuracy of these constants is not expected to be a limiting factor 
in any lattice implementation. This completes the renormalization of the theory at 
a level which will leave only 0(a) supersymmetry breaking effects. 

5. Conclusion 

Non-perturbative treatments of supersymmetric theories have been an elusive goal 
for decades. A recent surge in both interest and progress in this field has created a 
great amount of excitement in both the theoretical and lattice physics communities. 
However, these theories are only recently showing results in terms of the study of 
truly non-perturbative SUSY phenomena. Much of this difficulty has to do with the 
great technical challenges encountered in the attempt to formulate lattice theories 
that can reproduce SUSY in the IR. The direction of the field up until this point 
has been to try to build infrared supersymmetric behavior into the construction of 
the action. We have argued that, for three dimensional theories, it is feasible and 
straightforward instead to use the simplest possible action and do the fine tuning of 
its parameters necessary to obtain infrared supersymmetric behavior. This is possible 
analytically and does not prove as difficult as one might have feared. 

We have performed these tunings for a class of theories displaying N=2 super- 
symmetry in three dimensions and containing arbitrary numbers of matter multiplets 
transforming in the fundamental representation of the gauge group. The technique is 
robust in the sense that it relies only on theoretical principles, like Wilson's effective 
action formulation, and lattice implementations, like the Wilson action for fermions, 
that have been rigorously studied for decades. More generally, the entire lattice ac- 
tion is the most simple such construction with the appropriate IR limit, such that 
the many complications that can arise in such theories are suitably manageable or 
altogether absent. 

Though we have not done so, it should be very straightforward to extend our 
results to matter in general representations. This allows, for instance, the extension 
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of our results to SYM theories with higher supersymmetry, in which interesting non- 
perturbative physics is claimed to exist. For instance, the N=8 SYM theory in three 
dimensions has been conjectured by Seiberg to possess a non-trivial IR fixed point (an 
interacting conformal theory) ||37|| . The theory can be constructed as the N=2 SYM 
theory with 3 complex matter hyper-multiplets of the N=2 theory all transforming 
under the adjoint representation of the gauge group. A lattice has recently been 
constructed in |38| to study this theory by a very different technique. 
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A. Diagrams and Integrals 

In this appendix we enumerate precisely the integrals involved in the two-loop mass 
correction. All necessary diagrams are given in Fig. Heavy lines are fermions, 
light lines are scalars, dotted lines are ghosts and wavy lines are gauge bosons. Heavy 
dots represent the one-loop mass insertion counter-term, the one-loop diagrams eval- 
uated at zero external momentum. The dark cross in diagram (AX) is the measure 
counter-term. Diagrams that arise purely as artifacts of the lattice implementation 
are displayed in a dashed box. Diagrams that do not contribute in Landau gauge 
are not displayed. 

Since it makes structural sense to do so, we have grouped the expression for each 
two-loop diagram with the analogous mass insertion counter-term. The abbrevia- 
tions on the LHS correspond to the labels from Fig. |^. The notation is relatively 
straightforward: the first letter (F, S or A) labels the field which couple directly to 
the external scalar line (fermion, scalar or gauge boson respectively), the second, 
except in the case of the sunset diagrams and (AA6), labels the field coupling sec- 
ondarily and the number shows whether that secondary coupling is via a 3-point 
or 4-point vertex (replaced by C when the secondary line "crosses" the diagram). 
X represents the measure counter-term, a pure artifact of the lattice construction 
that comes about because the link variable integration measure in the path integral 



has non-trivial dependence on the gauge fields |3q| . G is the Fadeev- Popov ghost, 
which arises from gauge fixing in direct analogy to the continuum except that the 
ghost-gauge interaction contains an infinite number of new irrelevant vertices. Of 
these, only the GGAA vertex will contribute to continuum correlation functions at 
the perturbative order we require. Its contribution looks like a non-renormalizable 
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(SSUN) (SS4) (SF3) (FS3) (FSC) 



2-loop diagrams with fermions and scalars 




(ASUN) (SA4) (SA3) (FAC) (FA3) 




(AG4) (AF4) Pure Lattice Artifacts 



2-loop diagrams with scalars, fermions and gauge bosons 



Figure 2: Diagrams needed for the renormalization of SQCD. 

mass divergence which is essential to ensure that these divergences cancel in the 
correlation functions, as is guaranteed by gauge invariance. 

The integrands involving only scalars and fermions are (not including (SS4) since 
its contribution is exactly canceled by its analogous mass counter-term insertion) 

(SSUN) = A k A p A k+p 
(SF3) - c.t. = A 2 k Tr (s p S k+p - ^ + S 2 k+p )) 
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(FS3) - c.t. = -A fc Tr (s 3 p (S k+p - S p ) + S 3 k+p (S p - S k+P )) 

(FSC) = A k Tr (S p Sl +p ) . (A.l) 

For the scalar-gauge and ghost-gauge sectors we need to define the gauge boson 
propagator on the lattice (Landau gauge is used throughout the calculations): 

in terms of which they are 

(SA3) = ^Af A p A fc+p (A p + A fc+P )(2^5) M (55), 
(&44) - at. = ^Af (A^cosOv) - 1) + A^(oos(fc+p) M - 1)) ^ 
(ASWV) = A fe Af Ag p cos(^) M cos(^^^ 
(AS3) = Af Af A p A fe+p (2^5)^(5+A?) Q ^ 
(AS4) = -Af Af (a p cos(p m ) + A k+p cos(k+p)^5^ a 5 U f3 

(AGS) = AfAf A p A fc+p (SS M cos(|)^ a cos(^) a ^ 

(AG4) = lAfAf (Ajp>~, + At P 5^ M 5+S,)^^ • (A.2) 

For the fermion-gauge sector it is convenient to define the momentum part of the 
^i\)A and ipipAA vertices as 

V p = ct m cos (| V - ^"Pm and ^ = W + r cos (|)m ( A - 3 ) 
respectively. In terms of these the integrands take the form 

(FAC) = Af Tr (s 2 p V 2 U S l +P K + k) 
(FAS) - c.t. = ±Af Tr (sl(y» p+k S k+p V» p+k - \[V k ^S k V k v + ^-^l) 

+s 3 k+p (v» p+k s p v» p+k - \Ks k v k » + 

(FAA) - c.t. = ^Af Tr (s p 3 (V^ - r) + S fc 3 +p (\4 +2p - r)) 5, u 
(AF3) = Af Af Tr S p ^ +fe S fe+p ^ +fc <W 

(^4) = i Af Af Tr + S k+ Xjt +2p ) 5,J vP . (A.4) 

Finally, for the pure gauge sector, we have 
(AX) = Af Af S^a 5 u p = X -At 
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(AA3) = Af Af AH p 

x |-^ A (2p+A;) M cos(^)„ + 5^(2^+^)^ cos(|) A + 5 tiv {p-k) x cos(^|^) M 

x|+(5 /37 (2p+A;) a cos(-) / 3 - 5 Q7 (2A;+p)^ cos(-) 7 - 5 Q/3 (p- fc ) 7 cos(— -) Q | .(A.5) 

The diagram labelled (AA6) is an artifact of the lattice discretization. Its con- 
tribution can be solved for exactly in terms of the complete elliptic integral of the 
first kind. Borrowing some now standard notation from |29| we define 



1 W2 

-4 / d 3 x 



Ti 



tt/2 



3.17591153562522. 



(A.6) 



It is 



(AA6) = - 



1 r2 



4 



t A t~t D t D + -t A t D t A t L 
3 3 



k.p 



AfA'"' 

k V 



Repeated Lorentz indices are summed. t A represents the generator of SU(AT C ) in 
either the defining (T^) or adjoint (F^ c ) representation depending on the transfor- 
mation properties of the external scalar in question. Appropriate external indices for 
the product of fs is implied. The group factor in square brackets reduces with the 
standard relations to 



C F C 



or 



\(c A f 



for external $ a j 
for external 6 . 



The momentum integral is 



A fc A, 



k.p 



(k-p) 



k 2 p 2 



(A.7) 



where k - p = 4 £V sin(/cj/2) sin(pj/2). The cross terms from the square are odd 
in the integration variables and thus integrate to zero in the Brillouin zone. The 
numerator of the second term can then be rewritten as an average, since the integral 
is unchanged by pi <-> p 2 <-> P3, and we see that the exact answer 



k.p 



A^A^ 

k p 



A k / A 



4 / £ 
3 V4^ 



Last, though certainly not least, we have the diagram (AA4). We will not repro- 
duce the 4-point gluon vertex rule here but it is written out in |]36| (with the changes 
mentioned in Sec. [|). The diagram contains two pieces. One involves the group 
structure /abe/cde, as in the continuum, and contributes to Cg% nd and C a g l . There 
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is also a pure lattice artifact piece with group structure Tr t a t b tCt d pi us p er _ 
mutations. The contribution of this piece is separately gauge invariant and infrared 
finite; it can be computed in a way similar to diagram (AA6) above. 

In terms of these we need to do 9 numerical integrals. These are composed 
into IR safe gauge invariant combinations as described in Sec. |j. The numerical 
factors on the diagrams are a combination of many different factors. Some care must 
be taken in determining extra factors of (—1) from diagrams where e commutes 
across odd numbers of gamma matrices and where odd numbers of the conjugate 
generator — T* appear. This occurs most commonly in diagrams with crossed lines. 
A notable exception is the —2(SSUN) contribution in the first line of Eq. [A.9| but 



this contribution, when lines representing the auxiliary fields are included, also takes 
the form of a crossed diagram. 

For the correction to the singlet scalar the appropriate IR safe gauge invariant 
combination is (we have changed notation so that analogous mass insertion coun- 
terterms are included implicitly by the diagram labels) 

= A(SF3) - Q(FS3) + (FSC) - 2(SA3) + (FAC) + 2(FA3) + (FAA) . (A.8) 

167T 

For the correction to the fundamental/anti-fundamental scalars they are 



°9l 



3{SF3) - 7(FS3) + A(FSC) - 2(SSUN) - (SA3) + (FA3) + \(FAA) , 



16tt 2 v ' y ' y ' v ' v ' y ' 2 

{SF3) - 4{FS3) + {SSUN) - (AF3) - (AFA) - (ASS) + 2{AS4) 



°92 



16vr 2 
16tt 2 

syfund 

16tt 2 



2(SF3) - 6(FS3) + 3(SSUN) - (SA3) + 2(ASUN) + 2(FA3) + (FAA) , 

(SF3) - 2(FS3) + (FSC) - (SSUN) - (SA3) + (FAC) 
--(ASUN) + 2(FA3) + (FAA) - ^(AS3) + (ASA) - (AF3) 
-(AFA) + (AG3) - 2(AGA) + ^(AA3) + -(AAA) + 2(AX) , (A.9) 
and for the adjoint scalar, 

= A(SF3) - 8(FS3) + 2(FSC) , 

167T 2 

(jadj 

= A(SF3) - §(FS3) - (FSC) + A(SSUN) - 2(SA3) 

167T 2 

+ (FAC) + 2(FA3) + (FAA) , 

C adj 1 

^ = -4(FS3) + -(FSC) - (SSUN) - -(FAC) 

lo7r 2 2 

— (AF3) - (AFA) - (AS3) + 2(ASA) , 
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= -2(FS3) - -(FSC) + -(ASUN) + -(FAC) + 2(FA3) + (FAA) 
167r 2 2 2 

--(AS3) + (ASA) - (AF3) - (AFA) + (AGS) - 2(AGA) 
+ l -(AA3)+ l -(AAA) + 2(AX). 

(A. 10) 

These combinations are observed to be perfectly IR safe through numerical analysis. 
Note that, as mentioned previously, the label (AAA) contains only that part of the 



4-point gluon vertex from |36[ that reproduces the appropriate expression in the 
continuum. The other symmetric piece, described in Sec. f|just after Eq. [Q| , has a 
unique group structure and is handled separately. 
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